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Abstract 

Ab initio calculation of the hyperfine, P-odd, and P, T-odd constants for the 
YbF molecule was performed with the help of the recently developed tech- 
nique, which allows to take into account correlations and polarization in the 
outercore region. The ground state ( 2 Tji/ 2 ) electronic wave function of the 
YbF molecule is found with the help of the Relativistic Effective Core Poten- 
tial method followed by the restoration of molecular four-component spinors in 
the core region of ytterbium in the framework of a non- variational procedure. 
Core polarization effects are included with the help of the atomic Many Body 
Perturbation Theory for Yb atom. For the isotropic hyperfine constant A, ac- 
curacy of our calculation is about 3% as compared to the experimental datum. 
The dipole constant A& (which is much smaller in magnitude), though better 
than in all previous calculations, is still underestimated by almost 23%. Being 
corrected within a semiempirical approach for a perturbation of 4/-shell in 
the core of Yb due to the bond making, this error is reduced to 8%. Our value 
for the effective electric field on the unpaired electron is 4.9 a.u.=2.5 x 10 10 
V cm -1 . 
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a. Introduction. A number of papers [l|-[5| were devoted to the calculations of the P, T- 
odd interaction constants in the ground state 2 £i/2 of YbF molecule. These calculations are 
necessary for the interpretation of the ongoing experimental search for the Electric Dipole 
Moment (EDM) of the electron d e Q. It is expected that in this experiment it will be 
possible to put a more stringent bound on the EDM of the electron. However, in order 
to link experimentally measured P, T-odd frequency shift with d e , one needs to know an 
effective electric field on the uncoupled electron, which is characterized by a constant Wd- 
So, reliable calculations of this quantity are essential. The cited above calculations predict 
the values of this constant in the interval (—0.62 -= — 1.5) x 10 25 Hz e _1 cm -1 . 

It is known, that parameter Wd is sensitive to the spin density in the vicinity of the 
heavy nucleus (see, for example, 0). The same, of course, can be said about magnetic 
hyperfine constants. So, comparison of the calculated constants of the hyperfine structure 
on the m Yb nucleus with the experiment provides an important test of the quality of the 
Wd calculation. In previous calculations hyperfine constants were significantly under- 
estimated. In [[J, it was concluded that the spin-correlation effects of the unpaired electron 
with the outermost core 5s- and 5p-shells of ytterbium should be taken into account in order 
to perform more accurate calculations of the hyperfine and P, T-odd constants. After that 
it was shown that 4 /-shell can also contribute to the dipole part of the spin density M. In 
unrestricted Dirac-Fock (DF) calculations by Parpia ||, the values for Wd are in a reason- 
able (mutually consistent) agreement with our previous Relativistic Effective Core Potential 
(RECP)-based Q and semiempirical |],[| calculations. Recently a new technique to account 
for the most important types of the core- valence correlations and the core polarization effects 
with the help of the atomic Many Body Perturbation Theory (MBPT) was developed and 
proved to be very efficient for the calculations of the hyperfine structure and P, T-odd effects 
for the BaF molecule Below we report the results of the application of this method to 
calculation of the YbF molecule. The deviation of the calculated from the experimental 
value is analyzed and the final Ad magnitude is corrected with the help of the semiempirical 
procedure [Qj. 

b. Spin-rotational Hamiltonian. For the 171 Yb isotope, which has nuclear spin / = 
|, the molecular spin- rotational degrees of freedom are described by the following spin- 
rotational Hamiltonian (see @): 

H SI = BN 2 + 7 SN - D e nE + SAI 

+ W A k A n x S ■ I + (W s k s + W d d e )Sn. (1) 

In this expression N is the rotational angular momentum, B is the rotational constant, S 
and I are the spins of the electron and the Yb nucleus, n is the unit vector directed along the 
molecular axis from Yb to F. The spin-doubling constant 7 characterizes the spin-rotational 
interaction. D e and E are the molecular dipole moment and the external electric field. 
The axial tensor A describes magnetic hyperfine structure on the Yb nucleus. It can be 
determined by two parameters: A = (A\\ + 2Aj_)/3 and Ad = (A\\ — AjJ/3. The smaller 
hyperfine structure associated with the 19 F nucleus is neglected. The last three terms in ([!]) 
account for the P- and P, T-odd effects. First of them describes interaction of the electron 
spin with the anapole moment of the nucleus k A ||. The second one corresponds to the 
scalar P, T-odd electron-nucleus interaction with the dimensionless constant k$. The third 
one describes interaction of the electron EDM d e with the internal molecular field E mo1 : 
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W d d e = 2( 2 Z 1/2 \H d \ 2 Z 1/2 ). (3) 

These equations show that the constant \W d characterizes an effective electric field on the 
unpaired electron. Similar to the W d , the other P- and P, T-odd constants Wi depend on 
the electron spin-density in a vicinity of the heavy nucleus and reliability of their calculation 
can be also tested by comparison of the calculated and experimental hyperfine constants. 

c. Calculation scheme. The Generalized RECP or GRECP (see |T(J and references) 
calculation of the ground state 2 S of YbF molecule was performed by analogy with [Q. 
The main difference of the present calculation is that the pseudospinors corresponding to 
the 5s- and 5p-shells were frozen from the calculation of the Yb 2+ ion with the help of 



the level-shift technique (which is also known as Huzinaga-type ECP, see e.g. [ll|). It was 
necessary to freeze these pseudospinors, because polarization of the corresponding shells 
was taken into account by means of the Effective Operator (EO) technique (see [[| for 
details). RASSCF |~2| calculations with 5284 configurations were performed for 11 electrons 
distributed in RAS-1=(2,0,0,0), RAS-2=(2, 1,1,0), and RAS-3=(6,4,4,2) subspaces. The 
discribed above procedure follows the same lines as the one used in our previous calculation 
of BaF molecule §. 

The molecular relativistic spinor for the unpaired electron was constructed from the 
molecular pseudoorbital (p^ 1 

^ = E^M + EcM' mi=0 + ---, (4) 

i i 

so that the atomic s- and p-pseudoorbitals of ytterbium in (|J) were replaced by the un- 
smoothed four- component DF spinors derived for the same atomic configurations which 
were used in generation of basis s- and p-pseudoorbitals. The MO LCAO coefficients were 
preserved after the RECP calculations. As the spin-orbit interaction for the unpaired elec- 
tron is small, the "spin-averaged" valence atomic p-pseudoorbital was replaced by the linear 
combination of the corresponding spinors with j = I ± 1/2 (see || for details). 

EOs for the magnetic hyperfine interaction, for the EDM interaction (|2]) and for the 
anapole moment interaction were constructed by means of the atomic MBPT. The EOs 
include two most important correlation corrections, which involve all the core electrons. The 
first correction corresponds to the Random Phase Approximation (RPA), and the second 
one corresponds to the substitution of the valence DF orbitals by the Brueckner orbitals. 
The latter are found by solving the one-particle equations with the self-energy operator o~(e) 
added to the Dirac-Fock operator h DF : 

(hj) F + a(e n )) 4> n = e n (f) n . (5) 

For the operator a(e) we used diagrammatic technique developed in |T3[. In our calculations 



of this operator we have neglected excitations from the shells Is ■ ■ - 3d. Opposite to that, 
in RPA equations it was important to include all the core electrons. Both RPA equations 
and Brueckner equations were solved for a finite four-component basis set in the V N ~~ 2 
approximation (which means that SCF corresponds to Yb 2+ ), and matrix elements of the 
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EOs were calculated. The basis set included DF orbitals for Is . . . 6s, 6p, 5c? shells. In 
addition 7 — 18s, 7 — 18p, 6 — 17d, 5 — 17 / and 5 — 14g orbitals were formed by analogy with 



the basis set N2 of Molecular orbitals were reexpanded in this basis set to find matrix 
elements of EOs for the molecular wave function. 

The main advantage of this method is that there is no need to extend the active (valence) 
space in order to include core electrons. It is of particular importance when one is interested 
in such molecular properties as hyperfine, P-odd, and P, T-odd constants. Corresponding 
electronic operators are singular at the nucleus. For this reason, it is necessary to include 
all the core electrons in RPA equations. Corresponding extension of the active space would 
be extremely expensive. Another approach which allows core electrons to contribute to the 
spin density is the unrestricted DF molecular calculation (where polarization is taken into 
account U). 

d. Results. Explicit expression for the parameter Wd is given in (|j) and (]3|). Other 
electronic matrix elements, which correspond to the parameters A, and Wi of opera- 
tor ([!]) can be found in . All the radial integrals and atomic four-component spinors were 
calculated for the finite nucleus 171 Yb in a model of uniformly charged ball. It is well known, 
that atomic matrix elements of operator (|2]) are proportional to Z 3 . The same scaling is 
applicable to the constant Ws, while the matrix elements which contribute to the constants 
W\ are proportional to Z 2 . As far as the nuclear charge of the fluorine is 8 times smaller 
than that of the ytterbium, we have neglected contributions to the Wi parameters from the 
vicinity of the fluorine nucleus. The additional argument which justifies this approximation 
is that the unpaired electron in the YbF molecule is mainly that of the Yb atom, and thus 
the spin density is localized near the Yb atom. At present we do not have RPA program for 
the scalar P, T-odd interaction, so we have focused here on the calculation of the constants 
W d and W A - 

Results for the parameters of the spin-rotational Hamiltonian are given in table |. Com- 
parison of the results of the GRECP/RASSCF calculation with the results of the present 
GRECP/RASSCF/EO calculation confirms our previous conclusion that correlations (in- 
cluding spin polarization) with the core give a significant contribution to the hyperfine 
structure and to the P-odd and P, T-odd constants. 

Our final value for the hyperfine structure constant A differs by less than 3% from the 



experimental value WM. That means that the symmetrical part of the spin density in the 



vicinity of the Yb nucleus in our calculation is probably rather good. However, it is also very 
important to reproduce the asymmetrical part of the spin density, which accounts for the 
dipole constant Our value for A& is in a slightly better agreement with the experiment 
than the value from M, but is still underestimated by almost 23%. About a half of this 
difference can be explained by the fact that in our molecular calculation 4 /-shell of the 
Yb atom was frozen. It was first pointed out by Khriplovich that in the YbF molecule 
excitations from the /-shell can be important. In particular, they can explain the small 
value of the spin-doubling constant 7 ||15|| . It was shown recently that contribution of the 
/-shell excitations to the spin density can give significant correction to the constant A& P| . 
Using equations (19) and (31) from the paper Q we obtain the following estimates for the 
/-shell excitation contributions to A and A&: 

5A « -3 MHz, 5A d « 15 MHz. (6) 
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Note that these corrections arise from the admixture to the molecular wave function of the 
configuration with the hole in the 4/ shell. The weight of this configuration was estimated 
in || to be of the order of 4%. Such an admixture is a purely molecular effect and is not 
accounted for by the EO technique. So, we can conclude that 5A d can be added to our value 
for Ad, resulting in « 94 MHz, that is in a much better agreement with the experimental 
value 102 MHz. 

It is important that similar contribution of the 4 /-shell excitations to the constant Wd 
is strongly suppressed. Indeed, operator (0) mixes /- and ci-waves. For Yb atom, 4d-shell 
is very deep and its mixture with the 4/-orbitals by the molecular field is very small, while 
5<i-shell is weekly bound and does not penetrate into the core region. Similar contribu- 
tions to other constants Wi are negligible due to the contact character of the corresponding 
interactions. 

One can see that the values of the Wd constant from the unrestricted DF calculation || , 
the most recent semiempirical calculation || and the present GRECP/RASSCF/EO cal- 
culation are in a very close agreement now. It is also important that the valence electron 
contribution to the Wd in [||] is in 7.4% agreement with the corresponding RECP-based 
calculation || (see table |). Another recent DF calculation Q gives the value, which is two 
times smaller. 
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TABLES 



TABLE I. Parameters of the spin-rotational Hamiltonian for YbF. 





A 


^4d 


w d 


w A 




Method 


(MHz) 


(MHz) 


(1 25 Ji2-) 


(Hz) 


(kHz) 


Exper. 14] 


7617 


102 








Semiemp. [pj 






-1.5 


730 


-48 


GRECP/SCF (| 


4932 


59 


-0.91 


484 


-33 


GRECP/RASSCF § 


4854 


60 


-0.91 


486 


-33 


Semiemp. [|| 






-1.26 




-43 


DHF g 


5918 


35 


-0.31 


163 


-11 


DHF+CP @ 


7865 


60 


-0.60 


310 


-21 


DHF (rescaled) § 






-0.62 


326 


-22 


DF (unpaired elect- 












ron contribution) || 






-0.962 






Unrestricted DF g 






-1.203 




-22 


GRECP/RASSCF/EO 












(this work) 


7842 


79 


-1.206 


634 




GRECP/RASSCF/EO 












(with 4/-hole correction) 


7839 


94 


-1.206 


634 
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